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Abstract. An improved microscopic diffusion in stars is presented considering in detail the partly ionized stages of 
metals. Besides, the influence of degenerate electron-gas and of the contribution of radiation to the total pr essure 



has been ac counted for. The solution of the diffusion equations is then performed following the scheme of Thoul 



■it al. (1994). By defining one mean charged ion per element very few modifications are necessary to solve the 
improved diffusion scheme. (A portable FORTRAN routine is provided.) The change in the sound-speed profile 
of a solar model obtained with the new diffusion description is at most about 25% at r = 0.6 -Rq. The biggest 
effect on low-mass stars is expected near the turn-off, where the convective envelope is shallowest. However, only 
a difference of at most 40 K in the effective temperature could be observed when assuming either fully or partly 
ionized metals in the diffusion equation. Nevertheless, the surface metal distribution is strongly altered. 

Key words. Sun: interior - Stars: evolution - Stars: abundances - Hertzsprung-Russel (HR) diagram 



1. Introduction 



ct al., 2000] ). With the improved models globular clusters 



In the last few years the precise measurement of solar os- 
cillations has challenged theorists to compute solar mod- 
els of gradually higher accuracy. This demanded an im- 
prove ment of the existing input physics like equation of 



3n*e- 

state ( plogors ot al., 19 4Mj ), opacities ([glorias fc Rogers 



have been found to be about 1 Gyr younger than deter 
mined previously by models without diffusion. 

However, recent observa tions of the surface iron abun- 
dance of near tu rn-off stars ( Ramirez et al., 2001 ; Ramirez 



|1996| ) and nuclear reaction rates ( Adelbcrger ct al., 1998| ) 
In addition, the formerly neglected process of microscopic 
diffusion has been found to improve the agreement of so- 



fc Cohen, 2002] ) suggests that diffusion is much less efR- 



lar models with hclioscismic data considerably (Bahcall 
|fc Pinsonncault, 1992 ). Nevertheless, there exists a sig- 



nificant discrepancy in the sound speed just below the 
convec tive between the theoretical predicted and the ho - 



vecp 
sisiiii 



lioseis mic determined value, the reason of which is still 



cient in metal-poor stars with thin convective envelopes 
than theoretically predicted. Even worse, the models pre- 
dict an almost total depletion of heavy ele ment s in the sur- 
face of such stars at the turn off (see Sect. 3.2). Including 
the here neglected effect of radiative forces, the opposite 
effect can be obtained for some elements, e.g., Fe may 
then be enh anced considerably at the surface ( Richardj 
et al., 2002). Thus, radiative levitation is important in 



unknown. 

Unfortunately, the expression "diffusion" is not always 
used in the literature to cover the same physical pro- 
cesses). In order to avoid confusion, "diffusion " is defined 

Bahcall et al. (1998| ) or |5alaris 



in this work like, 



e.g. 



et al. (2000), to describe abundance changes due to pres- 



sure, temperature, and concentration gradients neglecting 
the effects of radiative forces. 

Beside solar models diffusion is now also implemented 
in many stellar calculations, e.g., for the computation of 



globular cluster isochroncs (Castellani ct al., 1997; Balaris 
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certain evolutionary phases, and should be included in fu- 
ture models. 

Nevertheless, by assuming an additional mixing pro- 
cess below the convective envelope it would be possible 
to inhibit any diffusion process. Such a process is also 
favoured to reduce the discrepancy in the sound speed just 
below the solar convective zone ( [Richard et ah, 1996 ). 

Various mechanisms have been proposed to cause ad- 
ditional mixing. Before, however, being able to determine 
the extent of additional mixing processes or other non- 
standard physics, diffusion should be followed accurately. 
Therefore, the validity of physical assumptions used to 
compute the diffusion efficiency is investigated in this 
work. A common assumption in the calculation of diffu- 
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sion constants is the complete ionization of all elements 
Basically in all stellar model computations including diffu- 
sion (e.g. Castcllani ct al., 1997 ; [galaris et al., 2000 ; Weiss 



& Schl jattl, 2000| ; |Chaboyer et al., 200l| ) this approxima- 
tion is made. Exception are the models of Richer et al 
(2000|)]or Richard et al. (2002) which account for partial 



ionization and radiative levitation in the diffusion treat- 
ment. 

Since the cross section for the main scattering process 
of ions in the stellar plasma, the Coulomb scattering, is 
proportional to the square of the ion charge, deviations 
from complete ionization may have an important influ- 
ence on the diffusivity of the elements. But the diffusion 
constant of a specific element is not simply direct propor- 
tional to the cross section, because diffusion has to obey 
the laws of mass and charge conservation. Thus, the cross 
section of each element has to be considered in relation to 
the ones of all other elements. An exact knowledge of the 
ionization stage of each element is therefore necessary. 

In order to obtain more accurate microscopic diffu- 
sion constants the assumption of fully ionized metals is 
dropped in this work. Instead, the ionization stage of each 
metal is considered in detail, where the ionization degrees 
of each element are determined by using an up-to-date 
EOS of |Irwin (200"o| ). 

In the next section the implementation of partly 
ionized elements into the solution of Burgers' equa- 
tions ( Burgers, 1969 ) is described. The changes in the 
solar sound-speed profile and in the evolution of metal- 
poor stars using this improved description are discussed 
afterwards. 



2. Partly ionized elements in Burgers' equation 

The method to calculate diffusion constants as described 



by |Thoulet al. (1994[ ), denoted TBL hereafter, is followed, 
where the basic concepts and quantities can be found. The 
aim is to improve this method in particular by includ- 
ing partly ionized atoms. In addition, a possible non-ideal 
electron gas and the contribution of radiation pressure to 
the hydrostatic-equilibrium equation are taken into ac- 
countQ, the effects of which are neither considered by TBL. 
Besides, in this section nuclear burning is neglected; it can 
simply be added afterwards^ 

2.1. Basic equations 

An overall neutral plasma consisting of S species is as- 
sumed, where a species is defined to be an element in one 
particular ionization stage. Electrons are supposed to be 
species No. S, and are usually not mentioned explicitly. 



1 A FORTRAN77 routine containing all the non-standard 
effects in the computation of the diffusion constants is publicly 
accessible under http://www.astro.livjm.ac.uk/~hs. 

2 In the models presented nuclear burning is, of course, fully 
included. 



Burgers' equation for mass conservation can then be writ- 
ten as (cf. Eq. (11) in TBL) 



dn s 
~~dt 



r 2 dr 



(r 2 u!,n s 



(1) 



where r is the radial distance from the centre, and n s 
is the number density of ion s. The diffusion velocity 
u) s has to be determined by solving the momentum- and 
energy-conservation equations (cf. Eqs. (18.1) and (18.2) 
of Burgers, 196£, neglecting magnetic fields) 



dP s p s dP g[ 



dr 
and 



+ 



P 



dr 



p es E = /i(K, m,u},r h ) 



5 p dlnT 



(2) 



(3) 



where P s , p s {— n s m s ), and p es (= n s q s ) are the partial 
pressure, mass density, and charge density, respectively, 
for ion s with m s being its mass and q s its charge; E is the 
electric field strength. The vectors m and u> contain the 
mass and diffusion velocity for each ion s. For a detailed 
definition of the friction-coefficient tensor K, the residual 
heat-flow vector r^, and the functions /i and /2 the reader 
is referred to TBL. 

The total pressure P is composed of radiation (P ra d) 
and gas pressure (-P gas ), thus 



P = P 



rad 



rad 



where the standard definition f3 — P &as /P has been em- 
ployed. Since the departure from local charge neutrality 
is very small (cf. TBL), p e E is close to zero, and thus the 
usual hydrostatic-equilibrium equation 



dP 

dr 



-gp 



(4) 



can be used, with g being the gravitational acceleration. 
Hence, Eq. (||) becomes 

^pL+/3g Ps -Bl^- Pes E = f 1 (K,rn,u>,r h ) (5) 
d?' p dr 

By two additional constraints, current neutrality, 
^q s n s uj s ^0, (6) 
and local mass conservation, 



m s n s u> s = 0, 



(7) 



Eqs. (^) and (|^) are completed to form a closed set of 
25 + 2 equations for the unknown quantities oj, r^, g and 
E. 

For a large number of species with varying abundances 
these equations are preferentially solved numerically. For 
this purpose, they are reformulated to be more suited for 
stellar evolution calculations. 
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The partial pressure P s for atom s is given by the ideal 
gas equation 



P s = n s k B T. 



(8) 



The electron pressure in stellar envelopes is well described 
by this formula, too. However, in the deep radiative re- 
gions of low-mass stars, and in particular in white dwarfs, 
the electron gas is degenerate. Therefore, an "effective" 
electron density is introduced 



Pe 

fc R T 



(9) 



such that the total pressure can be written 

where n$ — n e . Furthermore the concentration C s is de- 
fined by 



C s — — i 

Tip 



(10) 



analogously to TBL. 

With C — ^Cs one can now write 



(3P 



P, 



9l 
c 



and hence the first term of Eq. 



becomes 



dr 



(3P = CL / dln/3P + dln(C s /C) 
dr C \ dr dr 



(11) 



(12) 



Taking into consideration that in TBL the pressure P has 
been defined to be solely the gas pressure, this term seems 
to agree with TBL. However, here C s has been defined via 
Eq. (£o|) instead of TBL's Eq. (22) (C* s = n s /n e ), which 
accounts for a possible degeneracy of the electron gas. 

Demanding charge neutrality, J2 s =£e = — 9e^e, one 
gets rid of the dependence on one species, e.g. He ++ . Thus, 



c 



He++CHe+ + 



E <> c >)> 

s^Hc++,e 



(13) 



where £ s is the ionization stage of ion s ((^ s = for neutral 
atoms) and 



Up 



Similarly, the l.h.s. of Eq. (0) becomes 



5 dlnT 
2 s ^r~ 



5C. dlnT 
2 C P dr ' 



(14) 



(15) 



Setting Eqs. @-© into Eqs. d)-(@) one finds, after 
some algebraic manipulations, that the diffusion velocity 
can be written in general as (see Appendix |X|) 



T 5/2 



-6, 



(16) 




Fig. 1. a) The contributions of the six terms in Eq. fllTj ) to 
the diffusion velocity of hydrogen in a 1.1 Mq star near the 
TO with initial [Fe/H] = -1.3: £f (solid line), £j (short- 
dashed), £f (dash-dotted), 10 x £f (dash-dot-dot-dotted), 
10 x £| (dotted), and (long-dashed). The thin solid line 
represents b) The hydrogen mass fraction in this 

model. 



where 



Ai 



,dlnP 
dr 



■ dlnT 



+A^ +A c 
dr 



dr 

ding, 
dr 



E \. 

t^e,He++ 



E4, 



c dlnC f 
dr 

dlnCt 
dr 



(17) 



The diffusion coefficients A p , A T , A° , A^ ', and A? are 
obtained by the solution of Eqs. (||)-(0). For convenience, 
the six terms of Eq. @ are denoted £f , £f, £f, £f , C) 
and , respectively. Compared to the analogous equation 
in TBL (Eq. (40)), three additional terms appear here. 
They should not be considered as new types of diffusion 
but rather as corrections to pressure (term 4) and concen- 
tration diffusion (terms 5 and 6). 

In particular, the fourth term in Eq. ( |T7j ) (£^) cor- 
rects for the appearance of a radiation pressure in the 
hydrostatic-equilibrium equation. The fifth term (£f) 
takes into account that for degenerate electrons the ideal 
linear relation between pressure and density is not valid 
any more (Eq. (pT[)). Thus, this term is only important, 
when strong gradients in the electron degeneracy appear 
(e.g. in white dwarfs). 

In Fig. |l|a the contribution of each of the 6 terms in 
Eq. ([Tt]) to the diffusion velocity w of H in a 1.1 Mq 
star near the turn off (TO) are shown. Pressure (£f ) and 
temperature diffusion (£j) are clearly dominating in most 
parts of the stellar interior. Near the centre, where already 
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most of the H has been consumed, concentration diffusion 
becomes important. The contribution of remains 
very small throughout the star, while £J increases toward 
the centre, where the electron degeneracy is gradually in- 
creasing. 



The last term in Eq. ([L7|) which originates from 
the demand for charge neutrality (Eq. (|l3|)), needs a more 
thorough consideration. It would disappear, if each ion is 
treated as a separate species, because the ionization stage 
£ of each species, i.e., the number of electrons per ion, is 
constant. Charge neutrality is then, as in the case of com- 
plete ionization, automatically conserved by the definition 
of A c . But if, as shown below, a "mean" ion per element 
is defined which is in an average ionization stage £, the 
number of electrons per species is varying in the ionizing 
regions of a star. This would cause artificial concentration 
diffusion. Due to £f this effect is compensated, ensuring 
that only the total number density of a certain element is 
relevant for its concentration diffusion, independent of the 
respective ionization stages in two neighbouring layers. 

Note that using one mean charged ion per element can- 



not lead to ambipolar diffusion (see Babel fc Michaud 
199l["l}nd references therein). Hence the ionization state 
of each element is determined in this case solely by the 
local thermodynamic equilibrium and cannot be modified 
by diffusion. 

An example is provided in Fig. [l^,, where mean charged 
ions per element were used to compute the contributions 
to the diffusion velocity. In the layers near log 10 T w 4.4, 
just below the convective envelope (c.f. Fig. g), H is not 
completely ionized (£h ~ 0.995). This small deviation 
from complete ionization causes a large negative contri- 
bution of the concentration diffusion term (dash-dotted 
line) to wh, although the hydrogen abundance is constant 
in this region (Fig. |l]b). However, by the last term in 
Eq. (^) this artificial concentration diffusion is counter- 
balanced such that the sum of and £j remains zero in 
this area (thin solid line in Fig. [j]a). 

The diffusion of ions into regions, where they are out 
of thermodynamic equilibrium, yields an additional en- 
ergy sink (or source). However, the ionization of the most 
abundant elements, which are usually hydrogen and he- 
lium, leads in most relevant cases to the development of 
convective zones (Fig. |2|), where mixing occurs on much 
shorter timescales than diffusion. Thus, in the radiative 
layers only less abundant elements may be in differently 
ionized stages. The out-of-equilibrium distribution of their 
ions leads only to a small amount of energy compared to 
the total internal energy of the plasma, and is thus ne- 
glected here. 



2.2. Linear approximation for uj s 

Using a mean ionization stage for all elements, instead 
of treating each ion separately, is valid if the diffusion 




Fig. 2. a) The mean ionization degrees C, of H (solid line), 
He (short dashed), C (dash-dotted), O (dash-dot-dot- 
dotted), and Fe (long dashed) in the Sun, normalized to 



their respective charge number Z (obtained from [rwin's 
|2000| EOS). The envelope convective zone is located within 
the grey shaded area, b) Same as a) but for a 1.1 Mq star 
near the TO with [Fe/H] = -1.3. 



velocity u>\ is a linear function of the charge of ion i of 
element s. Then 



and therefore 



5>x = £ 



(18) 



(w s , + C s Suj s )n s (19) 



is obtained, where Q are the discrete ionization stages of 
s, and £s is the mean charge defined by ( s = Cs n l) / n s 
with n s = J2i n l- 

In order to verify to what extent the assumption of 
lu s being linear in £ is justified, the diffusion coefficients 
for a simple mixture of hydrogen, helium and iron with 
X = 0.76, Y = 0.23 and X(Fe) = 0.01 have been com- 
puted. Conditions similar to those at the base of the solar 
convective envelope were chosen, i.e., T — 2.24 x 10 6 K 
and p = 0.19 g/cm 3 . 

In Fig. ||a the symbols mark the pressure diffusion coef- 
ficients Ap of H, He, and Fe for different ionization levels 
of the latter, while H and He are assumed to be fully ion- 
ized. For Fe I the interactions can no longer be described 
by Coulomb scattering, but are due to the atomic polar- 
isability (see e.g. Michaud et al., 1978). It has been found 
that the diffusion coefficients of a neutral atom is about 
a factor of 100 higher than of its single charged ion (see 
ponzalcz et al., 1995 ). To obtain a similar value but us- 
ing the Coulomb scattering description an "effective" ion- 
ization degree of 0.1 for neutral atoms is adopted. This 
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2.0 



Fig. 3. a) Dependence of the pressure diffusion constant 
A P of H (A), He (O) and Fe (□) on the charge (() of 
the Fe-ion in a mixture of X = 0.76, Y = 0.23 and 
X(Fe) = 0.01 at T = 2.24 x 10 6 K and p = 0.19 g/cm 3 . 
Using a mixture of two different Fe-ions with charges £(i) 
and C(* + 1) the dotted line is obtained for the mean dif- 
fusion constants Ap ("exact value"), while the solid line 
represents the result employing the mean charge £ ( "linear 
approximation"). At ( « 1 A' p (He) changes its sign, b) 
The relative difference between Ap and A' P for H (dotted 
line), He (dashed) and Fe (solid). 



certainly only yields approximate values for the diffusion 
constants, and should be improved in future works. 

Obviously, Ap is not a simple linear function of £, 
but with increasing ionization the non-linearity becomes 
smaller. The dotted line in Fig. ||a shows the mean diffu- 
sion coefficient A p for a mixture consisting of two Fe-ions 
with number density rii and rii+i and discrete charges Q 
and Ci+ij respectively. If instead the mean ionization stage 
- = niC(«)+n i+1 c(«+i) j d a i,q fo u ows the solid line. 

In Fig. the relative error between the two values 
is plotted. Because of the large mean free path of Fel, 
Ap(He) has a different sign for neutral than for ionized 
Fe causing He to diffuse in the same direction as H. This 
behaviour can also be observed, if Fe is assumed to be fully 
ionized, and the ionization stage of He is varied (Fig. |J). 
Then Fe would diffuse in the same direction as H, when 
mostly He I is present in the plasma. 

But in real stars such "anomalous" diffusion would 
rarely occur, as the ionization of hydrogen and helium 
leads to the development of convective regions, where mix- 
ing is much faster than diffusion. In these regions metals 
are, at least partially, ionized, too (Fig. ||). Therefore, in 
the diffusion-dominated radiative layers helium and met- 



c 



Fig. 4. Same as Fig. || but for varying He-ionization de- 
grees. Note that at £ w 0.9 A'(Fe) changes the sign. 

als almost never are in their lowest ionization stages, and 
thus the difference between Ap and A' P remains small. 

In some cases, however, the He abundances at the 
outer layers are reduced considerably by diffusion such 
that no convection develops due to He ionization. Then 
the diffusion of each He-ion should be treated separately. 
But taking into consideration that the diffusion velocity of 
He I is not determined very accurately and the He abun- 
dance in that case is very low, using the mean ionization 
degree of He may still yield sufficiently accurate results. 

In summary, the diffusion of partly ionized elements 
in stars can in most relevant cases be determined by us- 
ing one "effective" ion for each species carrying the mean 
ionization degree £. With this simplification the compu- 
tation of the diffusion of partially ionized elements is not 
connected with larger computational time than using fully 
ionized particles and remains therefore feasible. 

3. Stellar and solar models 

For the calculation of the stellar models a descendant of 



the program by Kippcnhahn ct al. (1967 ) in the latest ver- 
sion described in Schlattl (2001) is used. In particular the 
opacities tables of [glesias & Rogers (1996) and the nu- 
clear reaction rates given by Adelberger et al. (1998) are 
employed. While previously the equation-of-state tables of 
Rogers ct al. (1996| ), known as OPAL- or ACTEX-EOS, 
has been implemented, now the analytic description of 



[rwin (200C) is applied. This step has become necessary, as 
the publicly accessible Rogers et al.'s ( 1996| ) tables do not 
provide quantities like the electron density or ionization 
degrees. [Irwin| 's ( p000[ ) equation of state (IRWIN-EOS) 
is a further development of the SIREFF-EOS, which it- 
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Table 1. Physics in solar models 



Model 


EOS 


Effects included 
non- ideal e~ P ra( j 


in diffusion 
partly ionized 


OP1 


OPAL96 








OP2 


OPAL01 








IR1 


IRWIN 








IR2 


IRWIN 


X 






IR3 


IRWIN 


X 


X 




IR4 


IRWIN 


X 


X 


X 



4xl(T 
3x10 



-3 



self was built from the EFF-EOS flEggleton et al., 1973| ). 
While EFF could be solved during the stellar evolution 
program ( "in-line" ) , IRWIN is pretty slow, so that I pro- 
duced spline tables on a pressure-temperature grid for dif- 
ferent hydrogen and metal fractions. The met al distribu- 
tion was chosen from Grevesse &: Noels (1993 ). 



The chemical evolution is determined by solving si- 
multaneously nuclear network and diffusion equation in 
a common scheme ( Schlattl, 1999 ), where convection is 
treated as a fast diffusive process. We follow the abun- 
dance changes of H, 3 He, 4 He, 12 C, 13 C, 14 N, 15 N, 16 0, 
17 0, 20 Ne, 24 Mg, and 56 Fe, which implies that 92% of the 
total metallic mass is considered. 

3.1. Solar models 

In my solar models the outer layers above an optical depth 
of 1000 have been taken from a 2D-hydrodynamical atmo- 
sphere model of Freytag et al. (1996| ) , which is determin- 
ing also the outer boundary condition of the interior solar 
model. These atmospheres extend down to regions, where 
the stratification is almost adiabatic, including the supera- 
diabatic layers just below the surface. By this means im- 
proved p-mode frequencies could be obtained like with for- 
merly applied ID-atmospheres (Schlattl et al., 1997), and 
the models become almost independent of the convection 



theory applied. For the latter the description of Canuto & 
Mazzitelli (1992j ) is used. As usual, the convection param 
eter acM> initial helium and metal content are adjusted 
to yield models with present solar luminosity, radius and 
metal fraction Z/X. In order to explore the effect of dif- 
ferent diffusion approaches these solar constraints are re- 
produced with an accuracy better than 10 -6 . A summary 
of all solar models computed in this work is provided in 
Table 

In a first step models using the diffusion description 
of TBL containing OPAL- or IRWIN-EOS have been 
produced. Recently, an updated version of the OPAL- 
EOS has become available (denoted OPAL01, hereafter). 
The main differences compared to the version of 1996 
(OPAL96) are the improvement of the activity expansion 
method for repulsive interactions (Rogers, 2001), the in- 
clusion of HcH + and He ++ in addition to H2, and H~ , 
and a relativistic treatment of the electrons. 

In Fig. [5|a the differences of the sound speed between 
the seismic model of Basu et al. (1997) and models with 
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Fig. 5. a) Difference in the run of sound speed between the 
seismic model of Basu et al. (1997) and solar models OP1 
(dashed line), OP2 (dash-dot-dot-dotted), IR1 (solid), and 
IR4 (dash-dotted), respectively, b) Sound speed of model 
IR2 (solid line), IR3 (dash dotted), and IR4 (dashed) com- 
pared to the reference model IR1. 



OPAL96-, OPAL01-, and IRWIN-EOS, respectively, are 
shown. Interestingly, model OP2 containing the updated 
OPAL01-EOS provides a worse agreement on average with 
the solar sound-speed profile than model OP1 with the 
older OPAL96-EOS. The main reason is that the adi- 
abatic inde x T\ decreases when t reating electrons rela- 
tivistically ( Bonanno et al., 2001 ). However, while the 
overall sound-speed profile has deteriorated, models with 
OPAL01-EOS provide a better agreement between the he- 
lioseismologically determined solar age and the mcteoritic 
age of the solar system ( Bonanno et al., 2002] ) 

Model IR1 was computed with the IRWIN-EOS, where 
electrons are also treated relativistically and where H2- 
and ITj~ -molecules are included. While in the central and 
convective region (r > 0.7 Rq), the sound speed of IR1 
and OP2 are very similar, between 0.25 < r/R Q < 0.6 
IR1 reproduces the solar-sound speed profile better than 
model OP2. It is beyond the scope of this work to analyze 
what causes the difference between OPAL01- and IRWIN- 
EOS. Anyway, with the latter the solar sound-speed profile 
can be reproduced very well, and it thus provides an at 
least equally good description of stellar plasmas as the 
OPAL01-EOS. 

Starting with model IR1 a series of solar models have 
been computed where subsequently non-standard effects 
have been added in the diffusion equation. The relevant 
physics included in each model is summarized in Table [l| 

Since the electron degeneracy in the Sun is very small 
(Fig. H), its implementation in the diffusion equation alters 
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Table 2. Characteristic quantities of solar models in this work. The indices i, s, c, and C.Z. denote initial, sur- 
face, core and bottom of convective envelope, respectively. acM is the parameter of |Canuto fc Mazzitclli ! s ( 1992| ) 
convection theory accounting for non-local effects, and is the total flux of neutrinos produced in the 8 B-decay 
and the 3 He(p,e + z^ e ) 4 He-reaction. In the last two columns the expected rates for the Gallium (Ga) and Chlorine (CI) 
experiments are provided, where 1 SNU (solar neutrino unit) corresponds to 1 event per 10 36 target atoms and seconds. 



Model 


Y t 


Zi 


oicm 


Y s 


Z s 


flo.z. 


Pc 


T a 




Ga 


Cl 




gem ° 


10 'K 


10 e cm- 2 s- 1 


SNU 


SNU 


OP1 


0.27478 


0.01992 


0.97803 


0.24429 


0.01807 


0.7131 


152.01 


1.5714 


5.090 


128.6 


7.624 


OP2 


0.27520 


0.01987 


1.0264 


0.24509 


0.01805 


0.7132 


152.02 


1.5711 
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Fig. 6. The solid line shows the ratio between the electron 
pressure assuming an ideal gas (n e kBT) and the actual 
electron pressure P e in model IR1, while the dashed line 
represents /3, the relative contribution of radiation to total 
pressure. 



the sound speed at most only about 0.01% in the centre 
(Fig ||d), where the degeneracy is highest. 

An even smaller change is obtained by taking into ac- 
count the radiation pressure which contributes at most 
about 0.1% at r = 0.65 Rq to the total pressure (Fig. |J). 
There a tiny change in the sound speed of about 10 -5 can 
be observed (Fig |5Jb), but overall it remains essentially 
invariant. 

With the inclusion of partly ionized elements in the dif- 
fusion treatment the sound speed has increased by almost 
0.02% at r = 0.6 R @ , but decreased by about 0.01% in 
the centre compared to IR3. Taking into account that the 
sound speed of IR3 (and equally IR1) has differed by about 
0.08% from the seismic model for 0.3 < r/R Q < 0.65, an 
improvement of about 25% could be achieved with model 
IR4. 

In the present stage, however, the discrepancies be- 
tween different input physics like equation of state are 
exceeding this value and thus the effect of partly ionized 
elements on diffusion appears to be negligible. But, in- 
cluding this effect reduces the uncertainty in solar models 



Table 3. Abundance variations of solar models of this 
work. The central abundance of 12 C was in all models 
0.7% of the initial one. 



Model 


4 He 


Surface (%) 

12 C 16 Q 


Z 


Center (%) 
4 He le O 


Z 


OP1 


-11.1 


-10.3 -10.1 


-10.1 


133.1 -2.2 


8.0 


OP2 


-10.9 


-10.1 -9.9 


-9.9 


132.7 -2.3 


7.9 


IR1 


-10.9 


-10.1 -9.9 


-9.9 


133.3 -2.3 


7.9 


IR4 


-10.9 


-10.1 -10.2 


-10.3 


133.2 -2.3 


7.9 



and thus helps to analyze the still existing shortcoming in 
model input physics. 

Turcottc ct al. (1998) has investigated the influence of 
a detailed treatment of the ionization stages in the dif- 
fusion equation, too. But in contrast to my results, they 
found a decrease of the sound speed when using partly 
instead of completely ionized elements. However, they fol- 
lowed a different strategy for computing solar models by 
using the same initial Z/X ratio in the models under 
scrutiny. Furthermore, they failed to compare directly two 
models where the ions are assumed to be either fully or 
partly ionized. Instead they compared their model, which 
contained a detailed treatment of the ionization stages, 
with one where the diffusion velocities of all metals heav- 
ier than O are assumed to be equal to the velocity of 
completely ionized Fe. Both differences might contribute 
to the discrepancies between our results. 

In Table | some characteristic quantities of the pre- 
sented solar models are summarized. Only very small dif- 
ferences can be observed in all quantities. With the up- 
dated OPAL01- and IRWIN-EOS the core temperature 
has been reduced a little compared to model OP1 with 
the OPAL96-EOS. However, with the inclusion of partly 
ionized metals in the diffusion equation the temperature 
increases weakly, yielding a somewhat higher value for IR4 
as for OP1, and thus slightly enhanced neutrino rates. 

The abundance variations of 4 He, 12 C, and 16 at the 
solar surface and in the centre compared to their initial 
values are summarized in Table |^. Considering the detailed 
ionization stages of the metals in the diffusion scheme yield 
a higher reduction of 16 O and all metals on average, but 
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Fig. 7. (a) The abundance changes at the surface between 
initial and solar model of sequence IR3 (dashed line) and 
IR4 (solid) . On the right hand side the change in the total 
metallicity Z is shown, (b) As (a) but in the centre. All 
elements lighter than O have been omitted, as their abun- 
dance is basically determined by the nuclear equilibrium 
of the CNO cycle. 
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Fig. 8. (a) The diffusion velocity of Fe in the Sun assum- 
ing fully (dashed line) resp. partly (solid) ionized elements 
and their relative difference (dashed line), (b) Like (a) but 
for O. In addition the relative difference in the diffusion 
velocity of C between part and complete ionization are 
shown, too (dash-dotted line). 



the absolute changes remain small (< 0.4%). The biggest 
change can be observed for iron, which is about 2.5% 
stronger depleted in model IR4 th an in IR3 (Fig. [^). Thi s 
is in agreement with the result of Turcotte et al. (1998| ), 



but the overall change in the metal fraction of 0.4% is 
about a factor of 2 less than found by them. The reason 
is a smaller change of the diffusion velocity in my com- 
putations by using partly instead of completely ionized 
metals. In particular, the velocity of oxygen, which is the 
most abundant element among the metals, is increased 
only by about 5% below the convective envelope (Fig. ^), 
while Turcotte ct al. (1998 ) found an about 10% higher 
value (see Fig. 13, therein). 

The cause of this difference is difficult to disentangle, 
but it might be that the ionization fractions obtained from 
the IRWIN-EOS do not agree with the ones of the OPAL- 



EOS. The latter has been used by Turcotte et al. (199q ). 
A further clue supporting this assessment is that the sur- 
face carbon abu ndance is almost una ltered in my models 
(Table ||), while Turcotte et al. (1998 ) found a depletion of 
8.5% instead of 8%. According to the IRWIN-EOS the ion- 
ization of C below the solar convective envelope is nearly 
complete (Fig. ||a), i.e., no change in the carbon deple- 
tion is expected, if the improved diffusion description is 
applied. 

Independent of whether complete or partial ionization 
is assumed in the diffusion equation, the somewhat dif- 
ferent diffusion velocities of each element cause a slight 
redistribution of metals, such that they are inconsistent 
with the metal composition assumed in the opacity ta- 



bles. Turcotte et al. (1998) has shown that the usage of 



about 0.1% compared to the Sun. They obtained this dif- 
ference in a model with elaborate treatment of the ion- 
ization stages. However, as mentioned above, the metal 
distribution in their model was modified stronger com- 
pared to the initial one than in model IR4. So, using 
monochromatic opacities probably had a bigger effect on 
their model, than it would have on model IR4. 

3.2. Low-mass stars 

In addition to the solar models also models for low-mass 
metal-poor stars have been computed. The aim is to exam- 
ine the influence of the improved diffusion scheme on their 
TO properties. Certainly the inclusion of radiative levita- 



tion can alter the models considerably ( Richard et al 
[2002 ), but it is nevertheless important to show how the 



more detailed treatment of ionization stages in the diffu- 
sion scheme modifies the stellar models. 

Basically, the same input physics has been used as 
for the solar models, where for the equation of state 



Irwin's (2000) description has been applied. Since the 2D- 



atmospheres implemented in the solar models do not suf- 
fice in temperature, gravity and abundances for metal- 
poor stars, the empirical T-r-relation of Krishna Swamy 
(1966) has been included. For the convection Bohm- 
Vitcnse's mixing-length theory was employed with 



monochromatic opacities may change the sound speed by 



a mixing- length parameter of 1.77, close to the calibrated 
value of a solar model containing the same physics. The 
metal distributions in the models were chosen to be ot- 
enhanced with [a/Fe] = 0.35. Consistent opacity tables 
have been produced for this mixture (A. Weiss, private 
communication) . 
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Table 4. Characteristic properties of low-mass stars at the TO with initial Z = 10~ 3 and Y = 0.24, whereas an 
a-enhanced mixture has been used with [a/Pe] = 0.35. AM c e . is the mass of the convective envelope. To show the 
changes in the abundances {x} = \og w (x/xi) has been defined, where the index i denotes the initial value. The first 
row for each mass represents the no-diffusion case, while the second and third show the values considering fully and 
partly ionized elements, respectively, in the diffusion equation. 
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Fig. 9. The evolution of 0.9, 1.0, and 1.1 M© stars with 
[Fe/H] = —1.6 (a-enhanced) treating all elements in the 
diffusion as partly (solid line) or fully ionized (dashed) . For 
comparison the evolutions without diffusion are shown, 
too (dash-dotted line). 

In Fig. |^ the evolution of 0.9, 1.0, and 1.1 Mq stars 
from the zero-age main sequence until the early-RGB are 
displayed in the HR-diagram. Their properties at the TO 
are summarized in Table |]. The inclusion of diffusion leads 
to smaller ages at the TO, which has also been found by 
various authors. The effect of decreasing the age is more 
pronounced for smaller masses which have longer main- 
sequence lifetimes. There, diffusion is capable of decreas- 
ing the core hydrogen content considerably, the amount 
of which basically determines the time spent on the main 
sequence. 

Strongly varying pressure and temperature gradients 



these two parameters can be observed in the evaluated 
models. Since the mass of the convective envelope of the 

1.0 Mq star is orders of magnitude smaller than in the 
0.9 Mq star, helium is much more depleted in the for- 
mer than in the latter case (Table g). However, although 
the extension of the convective envelope is smallest in the 

1.1 Mq star, more He remains on the surface than in the 
1.0 Mq case; the reason for this is the smaller lifetime of 
the 1.1 Mq star on the main sequence. 

The surface metal abundances follow the trend of He. 
However, their total depletion and distribution is strongly 
depending on the treatment of ions in the diffusion equa- 
tion. Assuming complete ionization leads in general to 
smaller diffusion velocities (cf. Fig||) and thus to a smaller 
decrease of the surface metal content. Already in this case 
the metal distribution may be altered up to a factor of 
2 (Table ||). Taking into account the different ionization 
stages some metals may even totally disappear from the 
surface. 

However, whether this also occurs in real stars is ques- 
tionable, as with decreasing abundance and increasing ef- 
fective temperature the influence of radiative levitation, 
which has been neglected here, becomes stronger. It would 
diminish the tendency of metals to settle in the stellar in- 
terior and might even lead to an enhancement of some 
elements on the surface (Richard et al., 2002). 

An additional source of uncertainty in these models 
is given by the strong deviation of the metal distribution 
from the initial one, being apparent in particular in models 
with detailed ionization treatment. With an assumed stan- 
dard a-enhanced mixture in the opacity tables, the opaci- 
ties obtained are not valid any more at the TO, which may 



lead to different effective temperatures. However, Salaris 



near tl Le surface of these stars cause an almost total de- 



pletion of helium and metals at the surface near the TO 
phase. Thereby the degree of the helium depletion depends 
on two parameters, the extension of the convective enve- 
lope, which provides the reservoir for inward migrating 
atoms, and the time the star needs to reach the TO, i.e., 
the time diffusion can operate. The competing effects of 



ct al. (1993) have demonstrated that the TO temperature 
of a M = 0.8 Mq star with a metallicity Z = 10~ 3 is de- 
termined basically by the opacities above about 12,000 K. 
Moreover, the opacities there are only sensitive to the con- 
tribution of C, N, O, and Ne to the metals, whereas the 
detailed distribution of these elements is only of minor 
importance. Artificially decreasing their abundances, i.e., 
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Fig. 10. a) The helium profile of a 0.9 (solid line), a 1.0 
(dash-dotted), and a 1.1 (dashed) M star at the TO in- 
cluding diffusion with elaborate treatment of the ioniza- 
tion stages. The thicker lines mark the envelope convective 
zones, b) The contribution of C, N, O, and Ne to the to- 
tal metallicity. The line styles agree with a), however, the 
thick lines correspond to models with complete ionization 
assumed. The dotted lines represent in both panels the 
initial values. 



underestimating their contribution, yields too high effec- 
tive temperatures. 

In the models of this work, where full ionization is 
assumed, the contribution of C, N, O, and Ne has been 
overestimated using standard a-enhanced opacity tables 
(Fig. 10 3), implying that the obtained effective tempera- 
ture are actually too small. In way of contrast, just the 
opposite behaviour can be observed in the models where 
the ionization stages are considered in detail. The result 
is that the actual difference in effective temperature be- 
tween models assuming either complete or partial ioniza- 
ti on would be smaller as quoted in Table ^. From the work 
of Salaris et al. (1993 ), I conclude that the reduction would 
be about 10 K. 

Thus, when the redistributed metal abundances are 
taken into account in the opacities, the effective tempera- 
ture of the 1.1 Mq star is increased only be about 20 K by 
treating ionization in the diffusion equation in detail. The 
effective temperature of the 0.9 Mq star remains almost 
unaltered. Probably, radiative levitation may slightly al- 
ter the effective temperature, but it is not expected that 
it alters the presented results by much more than 10 K. In 



fact, Richard et al. (2002) have found an increase of about 
10 K for a 0.8 M© star when including radiative levitation 
and using monochromatic opacities. For the age determi- 
nation of globular clusters the exact diffusion treatment is 
therefore only of minor importance, in particular, as the 



effect of the improved diffusion treatment is waning with 
decreasing TO mass, and thus increasing globular-cluster 
age. 

4. Summary 

Starting from Burgers' equation for a multicomponent 
fluid the effect of partly ionized metals has been con- 
sistently included in the diffusion equation following the 
procedure of TBL. As an additional feature, electron de- 
generacy and the contribution of radiation to the total 
pressure has been taken into account. It has been shown, 
that for most stars it is sufficient to use one ion per ele- 
ment carrying the mean charge instead of computing dif- 
fusion for each ion separately. By this means, the elabo- 
rate treatment of partly ionized metals remains feasible. 
With increasing stellar mass, and thus decreasing stellar 
convective envelope, it may become necessary to consider, 
at least, all He-ions separately. However, since these stars 
are usually hotter, radiative levitation becomes important, 
which then should be included, too. 

The effect of the improved diffusion treatment of partly 
ionized metals has been investigated in case of the Sun 
and the TO properties of low-mass metal-poor stars. In 
order to obtain the accurate ionization degrees of each 
metal species the equation of state of Irwin (2000) has 
been employed. The sound-speed profile in solar models 
obtained with this new EOS is in as good agreement with 
the helioseismic determined one as models with the up- 
dated OPAL01-EOS. 

An improvement of up to 25% between r — 0.25i?© 
and r = 0.65i?© could be achieved by including the effect 
of partly ionization into the diffusion equation. However, 
this value is presently still exceeded by the uncertainties 
in other input physics like the EOS. Nevertheless, the 
changes between models treating metals fully and partly 
ionized have been fou nd to be smaller than claimed by 
Turcotte et al. (1998 ). The origin of this discrepancy is 



not clear, but might be due to the use of different equa- 
tions of state. 

In low-mass metal-poor stars the improved diffusion 
description causes strong deviations in the surface metal 
distribution from the initial one. In addition, the deple- 
tion of metals is much stronger than in the case when full 
ionization is assumed. The consequent change in the ef- 
fective temperature is at most about 40 K for a 1.1 M© 
star diminishing with decreasing mass. 

The influence of radiative levitation on the diffusive be- 
haviour in particular of the metals has been neglected in 
this work. This may reduce the amount of depletion, and 
may further alter the metal distribution. Besides, since 
the effect of using partly instead of fully ionized metals in 
the diffusion equation is biggest for small convective en- 
velopes, a small amount of additional mixing below the 
convective boundary may reduce the effect of diffusion 
strongly. 

With the mass of the convective envelopes of metal- 
poor stars with M > 1 M© being very small 10~ 8 M©), 
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also the amount of mass loss may strongly modify the 
surface abundances. However, mass-loss mechanisms are 
only very poorly understood. Therefore, first of all, precise 
stellar models are needed to disentangle the extent of mass 
loss and to determine the influence of rotation. 
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Appendix A: 

The algebraic transformations of Eqs. (|^)-(0) are toilsome but 
mathematically not very demanding. Therefore, only the basic 
intermediate steps and the final results are described here. 

First, the r.h.s. of Eqs. (jE]) is converted, which consists ac- 
cording to Eq. of two terms, where the first one, d 1 ^ p , 
is almost in its final shape 

dln/3P dlnP din/? 



dr 



dr 



dr 



Using Eq. (|13[), the remaining one can be written as 



dln(C s /C) ( * C t \ din C t 



dr 



Chc 



E 



dr 

tj5He,e 

CtCt d In C t , >r^ CtCi dlnCt 



E 



ChcChc dr ^ ChcCho dr 



(A.l) 



where "He" denotes either of the helium ions (not the neutral 
He!), and 



C 



C s , for all ions 
C e , for electrons. 



The l.h.s. of the temperature equation (Eq. (g)) is given by 
Eq. ©. 

The equations for current neutrality (Eq. (n)) and local 
mass conservation (Eq. ([?])) become, using Eq. ([to|), 



\ " 



A s CsUs — 



and 

E 



(sC a UJ a = 0, 



(A.2) 



(A.3) 



where A a is the mass of ion s in atomic units. The system of 
equations (p|)-(|7|) can now be written in the form 



K 



dlnP 
dr 



+ <p- 



dln/3 
dr 



+ Vi 



dlnT 
dr 



+E 7ij 



dr 



§ -. 2iS+2 

+E^^ E ^» • 

3=1 J j=l 



where He is supposed to be species No. 2 and the electrons 
species No. 5. The quantities Ko, a, v, 7, A, and W are 
defined as in TBL with the following exceptions: 



1- Hi = 7? 
i = l,... 

2. An 



S (cf. Eq. (O)) 



for 



3. Ay = 



4. Wi 



AjCj for j = 1, . . . , 5; i = 25 + 1 (cf. Eq. ( A2)) 
GCj for j = 1, . . . , 5; i = 25 + 2 (cf . Eq. (fO|)) 

GCj for i = 25 + l;i = l,...,5 
-AjCj for i = 25 + 2; i = 1, . . . , 5 

K^heeE for j = 25+1 
K^I3n e m Q g for j = 25 + 2, 

/too l. 5 / 2 „-4 



2tt K B 



where mo is the atomic mass unit and Ko = 
Finally, the new coefficients tp and A are defined (see Eq. (A.l)) 
as 

Ci AiCi e . 

= -77 =- tor % = 1 , . . . , b 



and 



-T7-(«o-^)gS fori = l > ...,5-l 
for i = 5, . . . , 25 + 2. 



To account for the partly ionized metals in the particle in- 
teraction, TBL's Z s has to be substituted by £ s in the cross 
section a st (TBL's Eq. (8)), and in the Coulomb logarithm 
lnA st (TBL's Eq. (9)). 
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